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GPS RECEIVER USING DIFFERENTIAL CORRELATION 
FIELD OF THE INVENTION 

5 

This invention relates to the reception of a spread-spectrum signal and, in 
particular, to the reception of a spread-spectrum signal of a positioning 
system, the signal having a spreading code that identifies one of a number of 
signal sources. 

10 

BACKGROUND OF THE INVENTION 

The Global Positioning System (GPS) is a positioning system that provides 
accurate, worldwide, three-dimensional position Information to users with 

15 appropriate receiving equipment GPS includes a number of satellites, 
currently 24 satellites orbiting in six orbital planes with four satellites per 
plane, and a worldwide ground control and monitoring networic for monitoring 
the health and status of the satellites. Each of the satellites transmits a signal 
with a specific structure and comprising navigation data. GPS receiving 

20 equipment receives the signals from a subset of the number of satellites and 
determines the current position of the equipment from the received signals. 
The position of a GPS receiver may be computed when signals from at least 
four satellites are reliably received and measured. The transmitted signals 
have a direct sequence spread spectrum structure, also known as direct 

25 sequence code division multiple access (DS-CDIVIA) structure. DS-CDMA is 
one of two approaches to spread spectrum modulation for wireless 
communication of digital signals. In direct sequence spread spectrum, the 
stream of information to be transmitted is spread over a larger bandwidth 
than actually necessary for communications. This is achieved essentially by 

30 multiplying the data with a spreading sequence that has a much higher rate 
than the rate of the data stream. The rate of the spreading sequence divided 
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by the rate of the data is commonly known as spreading ratio or processing 
gain. The spreading code helps the signal resist interference and also 
enables the original data to be recovered if some spreading chips are 
unreliable. In GPS, a DS-CDMA structure is used with a 1023-'Chip spreading 
5 sequence and a period of 1ms. 

Superimposed on the spreading sequence is a 50 bits per second navigation 
data sequence that contains the information required for positioning. The 
positioning algorithm implemented by a GPS receiver is based on a 
10 determination of signal propagation delays by measuring the signal arrival 
times and by extracting satellite positions and transmission times from the 
navigation data sequence. The GPS receiver combines multiple distance 
measurements to triangulate the receiver position. 

15 The signal detection at the receiver is carried out by correlating the received 
signal sequence with a reference code. Each satellite spreads its message 
with its own spreading code, a so-called Gold code, such that different codes 
have a small cross-correlation with each other. In order to determine the 
signal delay, the receiver correlates the received signal with the Gold code 

20 for a given satellite and for different relative delays of the reference code. 
The correlation value has a well-defined peak when the delay of the correct 
reference code matches the actual delay of the received signal. At other 
delays, or for codes different from the Gold code of the satellite that has 
transmitted the received signal, the cross-correlation is very low, typically at 

25 least 24 dB lower than the peak value. When a peak is detected, the position 
of the peak is subsequently tracked by a delay-locked loop in order to 
determine the precise location of the loop. 

One of the most demanding tasks of a GPS receiver is the initial acquisition 
30 of the satellites. Without any prior knowledge, e.g. as to the approximate 
position, an autonomous GPS receiver may need to search through all 
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satellites, all possible code delays, and all possible frequency offsets in order 
to identiiy the Gold code and delay parameter of the received signal. Hence, 
the initial signal acquisition may require a multi-dimensional search, thereby 
rendering the initial acquisition a time-consuming task requiring many 
5 computational resources. Since the received signal level is very low, typically 
20-50dB below the receiver noise floor, the correlation samples typically 
need to be accumulated over rather long dwell times. A further complication 
is the presence of frequency offsets in the received signal caused by satellite 
Doppler effects and due to errors introduced by the local oscillators of the 
10 GPS receiver. 

It Is known to provide a priori knowledge to a GPS receiver in order to reduce 
the search space of the initial signal acquisition, e.g. by limiting the search to 
satellites expected to be visible over the horizon from the present expected 
15 location. For example, some GPS receivers are integrated into mobile 
terminals. Such devices may receive Information from a cellular 
telecommunications network, e.g. about the approximate location of the 
device, thereby assisting the Initial acquisition. 

20 However, it is a problem of the above prior art systems, that they require the 
presence of a priori information in order to speed up the acquisition. For 
example, not all cellular networks provide such information, and autonomous 
GPS receivers that are not integrated with a mobile terminal may not have 
access to such information. Furthermore, in areas without network coverage 

25 this information is not available either. 

Hence it is an object of the present invention to increase the efficiency of the 
initial data acquisition. 
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The above and other problems are solved by a method of receiving a spread- 
spectrum signal, the method comprising correlating a received spread- 
spectrum signal with a reference signal to detect the presence of one of a 
number of reference spreading codes; wherein the correlating further 
5 comprises performing at least one of the following steps resulting in a 
differentiated correlation signal: 

- differentiating the received spread-spectrum signal and the reference 
signal; and 

- differentiating the correlation signal; 

10 and wherein the differentiated correlation signal comprises a sequence of 
signal samples, each signal sample having a complex value. 

Consequently, the differentiated correlation signal is a complex signal, i.e. the 
samples of the differentiated signal have complex values including a complex 
15 modulus/norm and a phase/argument, thereby allowing a detection of both 
the time delay and a frequency error from the complex differentiated 
correlation signal. 

When the method further comprises detecting a frequency error of the 
20 received spread-spectrum signal from the differentiated correlation signal, the 
method provides an efficient estimation of the frequency error and, thus, an 
efficient and accurate signal acquisition. 

Furthermore, the method described herein provides a resource-efficient 
25 estimation of the frequency error that does not require additional search in 
the frequency space. 

Furthermore, the method described herein can be operated at a wide range 
of accumulation windows and demanding parameter settings. 

30 
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It has been realised by the Inventor that the computational complexity of the 
initial signal acquisition can be significantly reduced when the correlation is 
combined with a differentiation. 

5 It is a further advantage of the Invention that it provides a fast signal 
detection process. 

It Is a further advantage of the invention that it provides a detection method 
that Is robust against frequency offsets and insensitive to frequency 
1 0 instabilities of the received signal. 

According to a preferred embodiment of the invention, the method further 
comprises accumulating the differentiated correlation signal to obtain a 
correlation value. Consequently, the signal-to-noise ratio of the resulting 
15 correlation value Is Improved, thereby allowing reliable signal detection even 
in situations with a low signal-to-noise ratio (SNR) of the received signal. 
Reliable signal detection at low SNRs is achieved even without network 
assistance information, e.g. from a cellular telecommunications network. 

20 It is a further advantage of the invention that any phase rotation per chip 
does not accumulate during the accumulation of the correlation signal. 
Consequently, the correlation values are robust against frequency errors in 
the received signal. This has the advantage that the initial search space is 
reduced, since the search for any given spreading code may be limited to 

25 different code delays without having to take different frequency ofteets Into 
account. 

In embodiments of the Invention, the frequency error is determined from the 
correlation value. 

30 

In a preferred embodiment, the method further comprises determining a 
frequency compensation factor from the argument/phase of the correlation 
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value. Hence, the effect of the frequency enor on the subsequent processing 
of the signal can be reduced. 

In another preferred embodiment, the method further comprises determining 
5 a time delay of the received spread-spectrum signal from the complex 
modulus/norm of the correlation value. 

Hence an efficient estimate of the time delay and the frequency error are 
provided even in situations with large frequency errors. 

10 

In another preferred embodiment, the method further comprises de- 
spreading the received signal and extracting navigation data, e.g. the 
navigation message of a GPS signal, from the de-spread signal. 

15 In situations where the SNR of the received signal is low, the dwell time 
required for reliable signal detection is long. In the context of GPS, 
particularly low SNR are present when the GPS receiver is used in difficult 
reception conditions, e.g. indoors. The correlation signal may be 
accumulated coherently or non-coherentiy. In a preferred embodiment, the 

20 correlation signal is accumulated coherently by summing the complex 
product values, thereby further reducing the required processing time for a 
given SNR. in prior art systems, the presence of navigation data does not 
allow unlimited coherent accumulation due to bit transitions of the navigation 
data. It has been realised by the Inventors that a coherent accumulation of 

25 the differentiated correlation signal over long accumulation times provides a 
significant improvement of the reliability of the signal detection, in particular 
for signals with a low SNR. 

Hence, in a preferred embodiment, the accumulating comprises coherently 
30 accumulating the differentiated correlation signal- In a further preferred 
embodiment the received spread-spectrum signal comprises a digital 
information message encoded as bits, wherein bit transitions of the digital 
information message occur at predetermined transition time intervals; and 
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coherently accumulating comprises coherently accumulating the 
differentiated conrelation signal over a time interval that is longer than half the 
transition time Interval. Long coherent accumulation times larger than half the 
transition time interval between bit transitions of the navigation message are 
5 possible, because the accumulated differential congelation signal Is 
insensitive to the navigation message. Consequently, even longer 
accumulation times are possible, even larger than the transition time Interval, 
thereby Improving the reliability of signal detection, even under difficult 
reception conditions. 

10 

In a preferred embodiment, the method comprises 

- providing a plurality of reference signals modulated by said one of a 
number of reference spreading codes and delayed by respective 
relative code delays; 

15 - correlating the received spread-spectmm signal with the plurality of 
reference signals to obtain a corresponding plurality of differentiated 
correlation signals; 

- accumulating each of the plurality of differentiated correlation signals 
to obtain a corresponding plurality of correlation values for respective 

20 code delays; and 

- detecting a correlation peak in the plurality of correlation values to 
Identify a code delay of the received spread-spectrum signal. 

In another preferred embodiment, the method further comprises 
25 - correlating the received spread-spectrum signal with a plurality of 

reference signals, each modulated by a corresponding one of the 

number of reference spreading codes, to obtain a corresponding 

plurality of differentiated conrelation signals; 

- accumulating each of the plurality of differentiated correlation signals 
30 to obtain a corresponding plurality of correlation values for respective 

reference spreading codes; and 
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- detecting a correlation peak in the plurality of correlation values to 
identify a spreading code of the received spread-spectrum signal. 

Hence, an efficient method of detecting peaks in the congelation values for 
5 different code delays and/or different satellites is provided, thereby efficiently 
detecting the visible satellites and code delays/pseudoranges. 

In a preferred embodiment, the conrelating comprises differentiating the 
received signal and the reference signal to obtain a differentiated received 

10 signal and a differentiated reference signal; and determining the 
differentiated correlation signal from the differentiated received signal and the 
differentiated correlation signal. In another preferred embodiment, the 
correlating comprises determining a correlation signal from the received 
signal and the reference signal; and differentiating the correlation signal to 

15 obtain the differentiated correlation signal. Hence, the steps of differentiating 
and determining a correlation signal may be interchanged. 

Preferably, the received signal is provided as a sequence of received signal 
samples; and the correlating and differentiating steps are performed on the 

20 signal samples. In spread-spectrum communications systems the signal 

samples are also referred to as chips, i.e. the chips correspond to the bits of 
the spreading code. Since the chip rate Is much higher than the bit rate of the 
transmitted data, each data bit or data symbol Includes a large number of 
chips. In embodiments of the invention, the differentiation is performed on a 

25 single-chip time scale. The differentiation involves a delay-conjugate-multiply 
operation. In particular the differentiation comprises a delay of the signal 
samples by a predetermined number of chips, the predetermined number 
being less than the number of chips per data bit / data symbol. 

30 In a prefenned embodiment, the reference signal is modulated at a 
predetermined chip rate by a reference spreading code comprising a 
predetermined sequence of reference samples, also referred to as code 
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chips; the correlating comprises sampling the reference spread-spectrum 
signal resulting in a sequence of received signal samples; the reference 
signal comprises a sequence of reference samples; and the con^eiating 
comprises 

- correlating samples of the sequence of received signal samples with 
samples of the sequence of reference samples to obtain a sequence 
of correlation samples; and 

- accumulating samples of at least a sub-sequence of the sequence of . 
correlation samples to obtain a correlation value. 

In a prefenred embodiment, differentiating comprises multiplying a signal 
sample with the complex conjugate of the delayed signal sample, i.e. one of 
the preceding signal samples of a sequence of signal samples. In one 
embodiment, the delay corresponds to one spreading chip, i.e. the preceding 
signal sample is the immediately preceding signal sample, thereby 
simplifying the system design, since only one delay element is required. 

In other embodiments, the delay con-esponds to a plurality of spreading 
chips, thereby providing a more accurate estimate of the frequency error. 
20 Since the maximum detectable frequency error depends on the selected 
delay. In a preferred embodiment, the delay is selected as large as possible 
while still allowing the estimation of the largest expected frequency 
uncertainty. In some embodiments, the delay/differentiation time is 
determined as a function of the receiver quality. In particular, the frequency 
25 stability of the local oscillator (LO) of the receiver determines a maximum 
feasible differentiation time, and in embodiments of the invention, the delay 
may be selected to be close to or substantially equal to the maximum 
feasible differentiation time. For example, in a high perfonnance receiver In 
which the LO has a frequency stability of 4kHz, a suitable choice of delay 
30 time may be 30 chips, while a suitable choice for an entry-level receiver with 
a frequency stability of 40kHz may be 5 chips. The chip delay time in most 



10 
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embodiments is between 1 and 100 chips, more preferably between 2 and 50 
chips. 

In altemative embodiments, the differentiating may be performed in a 
5 different way. For example, in particular in situations where the frequency 
enror is relatively small, the differentiating may be perfomied by subtracting a 
preceding signal sample from the current signal sample. In this embodiment 
it is preferred that, after the subtraction, the absolute value of the subtracted 
signal is taken before the subsequent summation. 

10 

In another preferred embodiment, correlating samples of the sequence of 
received signal samples with samples of the sequence of reference samples 
further comprises 

- differentiating the sequence of received signal samples and the 
1 5 sequence of reference samples; 

- multiplying the differentiated sequence of received signal samples with 
the differentiated sequence of reference samples to obtain the 
sequence of correlation samples. 

20 In yet another preferred embodiment, correlating samples of the sequence of 
received signal samples with samples of the sequence of reference samples 
further comprises 

- multiplying the sequence of received signal samples with the 
sequence of reference samples to obtain a sequence of multiplied 

25 samples; and 

- differentiating the sequence of multiplied samples to obtain the 
sequence of correlation samples. 

In yet another preferred embodiment, spreading code is indicative of one of a 
30 number of signal sources, e.g. one of a number of space vehicles of a 
positioning system such as GPS. 
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The spreading code may be a pseudo-random-noise code. In a prefen^ed 
embodiment, the spreading code Is a Gold code. 

5 Further preferred embodiments are disclosed In the dependant claims. 

It is noted that the features of the method described above and in the 
following may be implemented in software and carried out on a data 
processing system or other processing means caused by the execution of 

10 program code means such as computer-executable instructions. Here and in 
the following, the term processing means comprises any circuit and/or device 
suitably adapted to perform the above functions. In particular, the above term 
comprises general- or special-purpose programmable microprocessors, 
Digital Signal Processors (DSP), Application Specific Integrated Circuits 

15 (ASIC), Programmable Logic Arrays (PLA), Field Programmable Gate Arrays 
(FPGA), special purpose electronic circuits, etc., or a combination thereof. 

For example, the program code means may be loaded in a memory, such as 
a RAM, from a storage medium or from another computer via a computer 
20 network. Alternatively, the described features may be implemented by 
hardwired circuitry Instead of software or in combination with software. 

The present invention can be implemented in different ways including the 
method described above and in the following, an arrangement, and a device, 
25 each yielding one or more of the benefits and advantages described in 
connection with the first-mentioned method, and each having one or more 
preferred embodiments corresponding to the preferred embodiments 
described in connection with the first-mentioned method. 

30 The invention further relates to an anrangement for receiving a spread- 
spectrum signal, the arrangement comprising correlation means for 
correlating the received spread-spectrum signal with a reference signal to 
detect the presence of one of a number of reference spreading codes; 
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wherein the correlation means Is adapted to generate a differentiated 
conrelatlon signal and comprises at least one of 

- means for differentiating the correlation signal; and 

- means for differentiating the received spread-spectrum signal and the 
5 reference signal; 

and wherein the differentiated correlation signal comprises a sequence of 
signal samples, each signal sample having a complex value. 

In a preferred embodiment, the arrangement further comprises means for 
10 detecting a frequency error of the received spread-spectrum signal from the 
differentiated correlation signal. 

In a preferred embodiment, the arrangement further comprises accumulator 
means for accumulating the differentiated correlation signal to obtain a 
1 5 correlation value. 

In a preferred embodiment, the arrangement further comprises means for 
determining a frequency compensation factor from the angle argument of the 
correlation value. 

20 

In a preferred embodiment, the accumulator means is adapted to coherently 
accumulate the differentiated correlation signal. 

In a preferred embodiment, the received spread-spectrum signal comprises a 
25 digital information message encoded as bits, wherein bit transitions of the 
digital information message occur at predetermined transition time intervals; 
and the accumulator means is adapted to coherently accumulate the 
differentiated correlation signal over a time interval that is longer than half the 
transition time Interval. 



30 
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In a preferred embodiment, the an^ngement further comprises means for 
providing a pluraiity of reference signals modulated by said one of a number 
of reference spreading codes and delayed by respective relative code delays; 
the conrelation means is adapted to conflate the received spread-spectrum 
5 signal with the pluraiity of reference signals and to generate a con^esponding 
plurality of differentiated con-elation signals; and the arrangement further 
comprises accumulator means for accumulating each of the plurality of 
differentiated correlation signals and for generating a conespondlng plurality 
of correlation values for respective code delays; and peak detection means 
10 for detecting a conrelation peak in the plurality of correlation values and for 
identifying a code delay of the received spread-spectrum signal. 

In a preferred embodiment the correlation means is adapted to correlate the 
received spread-spectrum signal with a plurality of reference signals, each 

1 5 modulated by a corresponding one of the number of reference spreading 
codes, and to generate a corresponding plurality of differentiated correlation 
signals; and the arrangement further comprises accumulator means for 
accumulating each of the plurality of differentiated correlation signals and for 
generating a corresponding plurality of correlation values for respective 

20 reference spreading codes; and peak detection means for detecting a 
correlation peak in the plurality of correlation values and for identifying a 
spreading code of the received spread-spectrum signal. 

In a preferred embodiment, the arrangement further comprises means for de- 
25 spreading the received spread-spectrum signal and for extracting infonmation 
data from the de-spread signal. 

In a preferred embodiment, the means for differentiating is adapted to 
differentiate the received spread-spectrum signal and the reference signal 
30 and to generate a differentiated received signal and a differentiated reference 
signal; and the conelation means is adapted to determine the differentiated 
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correlation signal from the differentiated received signal and the differentiated 
reference signal- 
in a preferred embodiment, the correlation means is adapted to determine a 
5 correlation signal from the received spread-spectrum signal and the 

reference signal; and the means for differentiating is adapted to differentiate 
the correlation signal to obtain the differentiated correlation signal. 

In a preferred embodiment, the means for differentiating a signal comprises a 
10 multiplier for multiplying a signal sample of the signal with the complex 
conjugate of a preceding signal sample. 

In a preferred embodiment, the reference signal is modulated at a 
predetermined chip rate by a reference spreading code comprising a 

1 5 predetermined sequence of code chips; wherein the arrangement comprises 
sampling means for sampling the reference spread-spectrum signal resulting 
in a sequence of received signal samples; wherein the reference signal 
comprises a sequence of reference samples; wherein the correlation means 
is adapted to correlate samples of the sequence of received signal samples 

20 with samples of the sequence of reference samples to obtain a sequence of 
correlation samples; and wherein the arrangement further comprises 
accumulator means for accumulating samples of at least a sub-sequence of 
the sequence of correlation samples to obtain a correlation value. 

25 In a preferred embodiment, the correlation means comprises 

- means for differentiating the sequence of received signal samples and 
the sequence of reference samples; 

- multiplier means for multiplying the differentiated sequence of received 
signal samples with the differentiated sequence of reference samples 

30 to obtain the sequence of correlation samples. 
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In a preferred embodiment, the correlation means comprises 

- multiplier means for multiplying the sequence of received signal 

samples with the sequence of reference samples to obtain a sequence 
of multiplied samples; and 
5 - means for differentiating the sequence of multiplied samples to obtain 
the sequence of correlation samples. 

In a preferred embodiment, the spreading code Is indicative of one of a 
number of signal sources, e.g. one of a number of space vehicles of a 
10 positioning system such as GPS. 

In a preferred embodiment, the spreading code is a pseudo-random-noise 
code. In a preferred embodiment, the spreading code is a Gold code. 

15 The invention further relates to a device comprising such an anrangement In 
one embodiment, the device Is a communications device, e.g. a mobile 
terminal with an integrated GPS receiver. The term communications device 
comprises any device comprising suitable circuitry for receiving and/or 
transmitting communications signals, e.g. radio communications signals, to 

20 facilitate data communication. Examples of such devices include portable 
radio communications equipment and other handheld or portable devices. 
The term portable radio communications equipment includes all equipment 
such as mobile telephones, pagers, communicators, i.e. electronic 
organisers, smart phones, personal digital assistants (PDAs), handheld 

25 computers, or the like. 

In other embodiments, the device Is a stand-alone GPS device. 



BRIEF DESCRIPTION OF THE DRAWINGS 

30 
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The above and other aspects of the invention will be apparent and elucidated 
from the embodiments described in the following with reference to the 
drawing in which: 

5 fig. 1 shows a schematic blocl< diagram of a GPS receiver; 

fig. 2 illustrates the search space for the initial signal acquisition; 

figs. 3a-b illustrate more detailed functional block diagrams of embodiments 
10 of the digital receiver block; 

figs. 4a-b schematically illustrate different embodiments of the ordering of 
differentiation and correlation blocks. 

15 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig. 1 shows a schematic block diagram of an embodiment of a digital GPS 
receiver. The GPS receiver comprises an antenna 101, e.g. a right-hand 
circulariy polarized (RHCP) antenna with substantially hemispherical gain 

20 coverage, for receiving the GPS radiofrequency (RF) signals of all satellites 
In view. The received RF signals are fed Into a pre-amplifier 102 which 
amplifies the received signals. There may be a passive bandpass pre-filter 
between the antenna and the pre-amplifier to minimize out-of-band RF 
interference. The amplified and optionally filtered RF signals generated by 

25 preamplifier 102 are then fed into a down-conversion and analog-to-digitial 
conversion block 103. In block 103, the signals are down-converted to an 
intermediate frequency (IF) using signal mixing frequencies from local 
oscillators (not explicitly shown). Block 103 further performs an analog-to- 
digital (A/D) conversion of the down-converted signal. Furthermore, block 

30 103 may comprise additional functions, such as an automatic gain control 
(AGO). The resulting digitized IF signal Is fed to a digital signal processing 
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block 104. Hence, block 104 receives a signal comprising signals from all 
GPS satellites in view. Typically, the signals are buried in the thermal noise 
at IF. At this stage, no demodulation has yet taken place, only signal 
conditioning and conversion to the digital IF. The digital signal processing 
5 block 104 comprises a digital receiver block, typically comprising a number of 
digital receiver channels. The digital receiver block may be implemented in 
one or more application specific integrated circuits (ASICs). The digital 
processing block further comprises a receiver processing block 106 that 
implements baseband functions such as loop discriminators and filters, data 

10 demodulation, phase-lock indicators, and/or the like. In one embodiment, the 
receiver processing block 106 is implemented as a suitably programmed 
microprocessor. The receiver processing block 106 further performs 
decision-making functions associated with controlling the signal processing 
functions of the digital receiver block. In particular, blocks 105 and 106 

15 perform the initial signal acquisition and the subsequent tracking of the 

received signals to obtain timing information for the subsequent calculation of 
the cunrent position. Block 106 further comprises the demodulation of the 
received signals to extract the navigation message. A more detailed 
description of examples of the above components of a GPS receiver may be 

20 found in e.g. "Understanding GPS, Principles and applications", E. D. Kaplan 
(ed.), Artech House, 1996. Embodiments of the signal acquisition functions 
performed by blocks 105 and 106 will be described in greater detail below. 

The GPS receiver further comprises a navigation processing block 107, e.g. 

25 a suitably programmed microprocessor or other processing means. The 
navigation processing block receives the extracted navigation data and the 
timing infomiation from block 104 and calculates a position based on the 
received data using algorithms known in the art of GPS receivers (see e.g. 
ibid.). The calculated position data is finally presented to the user via user 

30 interface 108, e.g. a display. It Is understood that, alternatively or additionally, 
the position infomiation may be used as an input to further data processing 
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functions. Hence, alternatively or additionally to user interface block 108, 
other interface functions may be provided. 

It is noted that some or all of the above digital signal processing and/or data 
5 processing functions may be combined in the same circuitry. For example, in 
some embodiments a single high-speed microprocessor supports the 
receiver, and implements the navigation and user interface functions. 

Fig. 2 illustrates the search space for the initial signal acquisition. In order to 
10 perform the navigation processing, the GPS receiver needs to identify the 
satellites from which data is received, and the receiver needs to determine 
the relative delay of the corresponding received spreading code. Once the 
GPS signals from sufficiently many satellites are acquired, they may 
subsequently be tracked, e.g. by a delay-locked loop. As mentioned above, 
15 In the context of GPS the spreading codes are so-called Gold codes (see e.g. 
R. Gold, "Optimal Binary Sequences for Spread Spectrum Multiplexing," 
IEEE transactions Info. Theory, Vol. 33, No. 3, October 1967, pp. 619-621). 

Without any prior knowledge, an autonomous GPS receiver needs to search 
20 through ail satellites, all possible code delays, and all possible frequency 

offsets in order to identify the Gold code and delay parameter of the received 
signal. Hence, as shown in fig. 2, the initial search space may be illustrated 
as a three-dimensional search space where the search is performed over: 

- Multiple satellites: In the following the satellites will be labelled n = 

25 1,"sNv, where Nv is the number of visible satellites. In GPS, Nv ^24. If 

no a priori information is available allowing to rule out some of the 
satellites as not currently visible, all 24 satellites may have to be 
searched. 

- The code delay ijj' = 1 of the reference code with respect to the 
30 received spreading code. In GPS, the spreading codes comprise 1023 

chips, i.e. A/f= 1023. 
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- The frequency offeet f*. In the following the frequency offset is 

assumed to be scanned in an interval [-fmax> fmaxl in steps of Af, i.e. fk - 

finax + k'Af. I( = 0 Nf= 2fmax/Af. 

5 In fig. 2, the search space is illustrated as a number of 2D search areas in a 
ij-fk plane spanned by the tj-skIs 203 and the f/raxis 202, where each plane 
con'esponds to a satellite. The 2D search areas are designated 201-1, 201-2, 
203-3,and ,201 -N, Each area is divided into a number of cells, each cell 
corresponding to a satellite n. a frequency offset fk, and a delay % i.e. each 
10 cell may be indexed by ^nj,/f>. 

In tiie following, the processing time required for searching all the (Nv Nt- 
cells of the above search space is estimated. For the purpose of this 
estimation it is assumed that for each cell, e.g. for cell 204, a proper 

1 5 spreading code is selected and delayed by an appropriate delay value ^. 
Furthemnore, it is assumed that a received signal is compensated for a 
frequency offset fk, i.e. de-rotated, and the frequency-compensated signal is 
correlated with the delayed reference spreading code. The correlation results 
for groups of chips are assumed to be accumulated coherently and the 

20 accumulated correlation value is compared to a threshold in order to identify 
a peak correlation value. The required accumulation time, or dwell time, for 
each cell is selected to be Td^ = G Tc, where Tc is the chip period, i.e. the 
duration of a spreading chip, and where G Is a target processing gain that 
depends on the received chip SNR ^and on the desired SNR of the resulting 

25 decision variable. For the purpose of the current estimate, the target 

processing gain is assumed to be lOdB. In that case ^CVychips have to be 
accumulated coherently, in order to rise lOdB above the noise floor for 
making a decision. 
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Furthermore, the resolution 2lf of the frequency grid 202 also depends on the 
dwell time, laecause better frequency stability is required for longer 
accumulation times. For coherent accumulation, the maximum possible dwell 
time for a given frequency resolution may be approximated as Taw = 2/(3Af). 

5 

Accordingly, under the above assumptions, the total processing time Tr of 
the initial search is 

Tr = (Ny ' Nt ' Nf) Taw = 3 Nv M • W (10 Tc/ yf. 

10 

Even though, the search process may be heavily parallelised, the above 
linearised processing time Tr serves as a good measure for the complexity of 
the Initial data acquisition method. 

15 In particular, it is noted that the processing time for a complete search of ail 
the cells of the above search space is proportional to the size f^ax of the 
frequency search space, thereby particularly extending the processing time 
for difficult reception conditions. When the frequency uncertainty W is large, 
and when the SNR is low, the required dwell time is large. 

20 

It is further noted that in prior art systems the coherent accumulation times 
cannot be arbitrarily increased. In GPS signals, the navigation messages are 
encoded at a rate of 50bps, i.e. there is a bit transition every 20ms. In an 
unassisted situation, i.e. without a priori ioiowledge about the bit boundaries 

25 and the message contents, the maximum length for coherent accumulation 
has to be smaller than 20 ms, e.g. no longer than 10ms. However, when e.g. 
a 100ms coherent accumulation is replaced by a coherent accumulation of 
10ms blocl<s and subsequent non-coherent accumulation of the accumulated 
bocks, 10^ = 100 blocks need to be accumulated in order to achieve the 

30 same SNR. Consequently, the processing time in prior art systems grows 
increasingly as the SNR degrades. 
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Additionally, the long accumulation times needed to work at low SNRs before 
dispreading require very narrow coherence bandwidths, i.e. a high stability of 
the local oscillators. Therefore, at low SNRs the obtainable stability of the 
5 local oscillators may become the limiting factor. 

With reference to figs 3a-b and 4a-b, a method will be described that 
improves the efficiency of the initial signal acquisition and, thus, the 
processing requirements. 

10 

Figs. 3a-b Illustrate more detailed functional block diagrams of embodiments 
of the digital receiver block. The digital receiver block 104 receives the 
digitized IF signal from the down-conversion and A/D block 103 shown In fig. 
1 . Hence, the received signal may be modelled as a sequence of chips 
15 according to the chip rate of the spreading code. The received signal at chip 
time k may be modelled as 

Hence, the received signal is modelled as a sum of contributions from each 
of the A/v visible satellites indexed /7=^,...,A/v,and as an additional receiver 
20 noise component Vk. Here, gn denotes the path loss from the n-th satellite, 

b^"^ denotes the navigation message bit from the n-Xh satellite at time k, i.e. 
fjin) _ j.^ ^ c^") denotes the chip value of the spreading code conresponding 
to satellite n at time k, and g)„ denotes the phase rotation per chip for the 
signal from satellite n. 

25 

Referring to fig. 3a, the received signal is fed into the digital receiver 
channel block 105 which determines accumulated correlation values z(n,d), 

d=0,. . .,1023, n = 1 Nv. for each of the satellites and for different code 

delays d. In the functional block diagram of fig. 3a, this is illustrated by the 
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plurality of correlation and accumulation blocks 310, each of which 
calculating the correlation values for a given satellite and for different code 
delays. Accordingly, each of the correlation and accumulation blocks 310 
comprises 1023 multiplier blocks 304 and 1023 summation blocks 305. It is 
5 understood, however, that the actual degree of parallellsation of the 
correlation and accumulation functions depends on the actual 
implementation. The digital channel receiver block further comprises a 
differentiation block 306 that differentiates the received signal and feeds 
the differentiated signal into the correlation and accumulation blocks 310. 

10 The differentiation block 306 comprises a delay element 301 for delaying the 
received signal by one chip, a conjugation block 302 for generating the 
complex conjugate of the delayed signal and a multiplier 303 for multiplying 
the received signal with the complex conjugate of the delayed signal. Each of 
the correlation and accumulation blocks 310 comprises multiplier blocks 304 

15 for correlating the differentiated received signal with corresponding delayed 

and effective reference spreading codes s^"^^ = summation 

blocks 305 for accumulating the respective correlation signals. The effective 

reference spreading code s\^2^ for a given delay d and satellite n may be 

generated directly In a code generatton circuit (not explicitly shown) or 

20 generated by a differentiation block from the original spreading codes c^j^}^ , 

which in turn may be generated by a code generator, as illustrated in fig. 4a. 
Code generators for spreading codes are known in the art of GPS receivers. 
Each summation block accumulates the signal over a predetermined number 
of K chips, l-ience, the resulting correlation value z(n,cl) may be written as 



25 



/c=1 

The number K of chip values to accumulate depends on the assumed SNR of 
the received signal. In one embodiment, a signal acquisition will be 
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performed for a predetermined value of K. If the acquisition fails, K is 
incrBased, and the acquisition process Is repeated. The preferred numerical 
value of K depends on the signal-to-noise ratio. In one embodiment, K may 
be selected to be /C=1023. 

5 

In the following, the processing of the differentiation block 306 and the 
con-elation and accumulation block 310 will be described In greater detail. For 
the purpose of the following description, the index n will be disregarded, 
con^esponding to a case where A/v=1 , i.e. the case of a single satellite. This is 
1 0 merely done to ease the readability of the following description and is not 
Intended to limit the scope of the invention. A skilled person will Immediately 
appreciate that the following results are equally applicable to a multi-satellite 
situation where cross-correlations to unrelated codes are negligible, as is the 
case for the spreading codes used in current GPS implementations. 

15 

Hence, In the simplified case of A/v=1, the correlation value may be written as 

= ^bk-Abk-A-i^'k-ACk-A-^to^-d^^-d-^h^^ +V'k_^+VkV*k_^. 

k=^ 

Defining a new effective spreading code Sk according to Sk = CkCk-'t . the 
above expression may be rewritten as 

, K 

20 z(d) = eJ^ ^bk-Ab„_A-^Sk-ASk-d 

where w is the total noise term after differential accumulation. 

Consequently, the correlation operation Is effectively perfonnried with respect 
to the above new spreading code Sk. It Is noted tiiat Sk is a Gold code of the 
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same ortler and with the same autocorrelation properties as the original Gold 
code Ck. 

it has further been realised that the above con-eiatlon value is substantially 
5 insensitive to the navigation message bk, since the product bubk.i=1 except 
when the bit transitions occur. However, since the bit transitions occur only at 
a small fraction of the time. e.g. for of the order of 1 out of 1 0® chips, the enror 
introduced by the bit transitions is negligible. 

10 Hence, it Is an advantage of the invention that it provides a detemilnation of a 
decision variable that Is insensitive to the navigation message, thereby 
allowing long coherent accumulation times. This in tum has the advantage 
that an efficient signal acquisition is provided even in the case of difficult 
reception conditions. 

15 

It has further been realised that the phase rotation term merely contributes as 
an overall factor, I.e. the phase rotation is not cumulative. Hence, it is an 
advantage that the calculation is not sensitive to frequency offsets or stability. 

20 Consequently, the search space for the Initial signal acquisition is 
considerably reduced, since a search along the firequency axis Is not 
required. This is a particular advantage in situations with difficult reception 
conditions where the frequency uncertainty fmax is large. 

25 It Is a ftjrther advantage that long accumulation times may be achieved, 

thereby allowing operation at low SNR, e.g. during Indoor reception, without 
imposing high stability requirements on the local oscillators, thereby reducing 
the cost for the required components. 

30 In comparison, the direct correlation of the received signal with the spreading 
code according to 
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/c=1 1^=1 

depends on the actual values bk of the navigation message and the phase 
rotation e'''' provides an accumulative phase rotation over K chips that cannot 
be neglected unless the residual phase tp is sufficiently small. 

5 

The linearised complexity of the initial search in the system of fig. 3b is 

TD=N^NtT^Z^=NyNt^Ta. 

Y 

i.e. independent of the magnitude W of the frequency uncertainty. Here 
10 is the dwell time. 

it has been observed that the total operation count for the search can be 
reduced by a factor of 1 .5-2, depending on the frequency uncertainty. 

15 Still referring to fig. 3a, the conflation sums z(n,d) generated by conelation 
and accumulation bloclcs 310 according to eqn. (1) are fed into the 
processing blocl< 106. The processing blocic comprises a peak detection 
module 308 for detecting peal< delay values d'n for one or more satellites, e.g. 
by comparing the largest |z(>7,c/VI for each satellite n to the conrelation floor 

20 for c/i^'n. 

In the embodiment of fig. 3a, the code delay values d'n for one or more 
satellites are fed into the navigation processing block 107. In a network- 
based mode of operation, the code delays may be utilised as pseudoranges 
25 as an Input to the network-assisted determination of the current position, see 
e.g. Heikki Kaaranen et al., "UMTS Networi<s: Archtectures, Mobility and 
Services," John Wiley & Sons, 1'* edition, August 2001, ISBN 047148654X, 
p. 164. 



wo 2005/060118 



PCT/EP2004/014372 



26 



Fig- 3b shows a functional block diagram of an alternative embodiment 
similar to the embodiment of fig. 3a. However, this embodiment is well-suited 
for an autonomous GPS receiver or a GPS receiver operating In an 
5 autonomous mode. In this embodiment, the processing block 106 further 
comprises a frequency detection module 309 which detemriines the angle 

argument ^^''^ of the peak correlation sum according to 
^^'^)=Zz{n,d;,). 

10 In the embodiment of fig. 3b, the frequency detection module 309 generates 
a number of phase error estimtes, one for each of the visible satellites 
detected by the peak detection module 308. Furthermore, the processing 
block 106 outputs the corresponding delay values d'n, n=1,...,Na for the Na 
detected visible satellites. The digital signal processing block 104 of the 

15 embodiment of fig. 3b further comprises a number of correlation blocks 307, 
one for each of the Na detected visible satellites and each receiving a 
corresponding delay value and a con-esponding phase error estimate. Each 
of the conrelatlon blocks 307 determines the correlation of the received signal 
to the determined satellite code at the corresponding detemiined code delay 

20 according to 

k=^ 

Here denotes the m-th bit of the navigation message from the n-th 

satellite, and the correlation is performed over L=20-1023 chips, 
corresponding to the duration of one bit of the navigation message 
25 transmitted at a bit rate of 50 bits per second, and resulting in a frequency- 
compensated correlation. Accordingly, each correlation block 307 comprises 
a phase factor determination block 312 for determining a suitable phase 
factor corresponding to the detected phase error, and a multiplier 313 for 
multiplying the received signal with the determined phase factor. The 
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correlation block 307 further comprises blocl^ 31 1 which generates the 
delayed spreading code for the corresponding satellite n, delayed by the 
detemnined delay dn- The output of blocks 311 and 313 are multiplied by 
multiplier 314, and the result is fed into accumulator 315. The accumulator 

5 315 generates the above frequency-compensated correlation which is 
fed into the navigation processing block 107. 

The navigation processing block 107 comprises a number of navigation 
message recovery blocks 316, one for each visible satellite. Each navigation 
1 0 message recovery block 31 6 recovers the corresponding navigation 

message from the conrelatlons y{^^ . The data from the navigation messages 

from different satellites are then fed into location calculation block 317 that 
solves the corresponding location equations to determine the current position 
of the receiver. 

15 

In the embodiments of figs 3a-b, the differentiation of differentiation block 306 
is performed on the incoming signal. In altemative embodiments the 
differentiation is performed after multiplication of the received signal with the 
reference sequence, as will be illustrated with reference to figs. 4a-b. 

20 

Figs. 4a-b schematically illustrate different embodiments of the ordering of 
differentiation and correlation blocks. 

Fig. 4a illustrates the differentiation block 306 and one of the subsequent 
25 multipliers 304 and corresponding accumulators 305, as were described in 
connection with fig. 3a. The arrangement of fig. 4a further comprises a 
reference differentiation block 400 that generates the differentiated and 
delayed reference spreading code Sk-d from the corresponding delayed 
original spreading code c/^^. The delayed original spreading code c/^-tf is 
30 generated by a code generator circuit 410. The differentiation block 400 
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likewise comprises a delay block 401 for delaying the spreading code Ck-a by 
one chip to obtain ck^i, a conjugation block 402 for generating the complex 
conjugate c\^,i of the delayed code, and a multiplier 403 for multiplying the 
original code chip with the delayed code chip to obtain Sk^= Ck-<ic*k^i. It is 
6 noted that, for spreading codes comprising real numbers, the conjugation 
block may be omitted. 

Hence, In the embodiment of fig. 4a, the samples of the sequence of 
received signal samples are correlated with samples of the sequence of 
10 reference samples by 

- differentiating the sequence of received signal samples and the 
sequence of reference samples; and 

- multiplying the differentiated sequence of received signal samples with 
the differentiated sequence of reference samples to obtain the 

1 5 sequence of correlation samples. 

It is understood that the receiver may directly generate the effective codes Sk 
rather than first generating and subsequently differentiating the conventional 
codes c/c. 

20 

Fig. 4b illustrates an alternative arrangement in which the differentiation is 
performed after multiplication of the received signal with the reference 
sequence. In this embodiment, the received signal and the delayed reference 
spreading code are correlated by multiplier 407, and the correlated signal is 
25 fed into the differentiation block 408, comprising a delay block 404, a 

conjugation block 405, and a multiplier 406 as described in connection with 
the previous differentiation blocks. The differentiation block 408 generates a 
differentiated correlation signal which is fed into accumulator 305 that 
generates the correlation sum z(n,d) as described above. 
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Hence, in the embodiment of fig. 4b, the samples of the sequence of 
received signal samples are con^elated with samples of the sequence of 
reference samples by 

- multiplying the sequence of received signal samples with the 

5 sequence of reference samples to obtain a sequence of multiplied 

samples; and 

- differentiating the sequence of multiplied samples to obtain the 
sequence of correlation samples. 

10 The relation between the embodiments of fig. 4a and 4b may be expressed 
as follows 

(Ck-dTH) (Ck^itk-i)* = (Ck-dC*k-d-i) (rKr%i)> 

15 It is noted that the output of the correlators 304 and 406 and, thus, the 

accumulated correlation value z(n,d) are complex values, I.e. z(n,d) has the 

form z{n,d) = pe^^ , where p=|z(n,d)| is the complex modulus/norm of z and 

^ is the phase/argument of z. As described above, the frequency error is 
determined from the phase of the correlation value z and the time delay is 
20 determined from the complex modulus/norm p. 

In the above embodiments of a GPS receiver have been described that 
provide an efficient initial signal acquisition and reduced processing 
requirements. The GPS receiver described above implements a signal 

25 acquisition method that is not affected by the frequency offset and stability 
and is insensitive to the navigation message. In particular, the above method 
and arrangements for detecting correlation peaks in a received signal do not 
require a scan over all possible frequency offset values. The frequency 
uncertainty is removed via a differential operator and accumulation of 

30 differentiated correlated chip values. As a result practically unlimited coherent 
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accumulation may be performed without requiring strict piiase/frequency 
synchronisation of the receiver and without requiring knowledge of the 
navigation message, since the per-chip rotation remains negligible. The 
differential correlation sums are used to detect the correct code delays and 
5 the correlation sum corresponding to the true delay yields the frequency 
offset value. The receiver may use the detected peal^ position, i.e. the 
detected code delay, and the determined frequency offset value to perform 
conventional non-differential decon-eiation to recover the navigation message 
in an autonomous or a networlc-assisted mode. In a network-based mode, 
1 0 the code delay values may be reported as pseudoranges to the network 
directly. 

Although preferred embodiments of the present invention have been 
described and shown, the invention is not restricted to them, but may also be 
1 5 embodied in other ways within the scope of the subject matter defined In the 
following claims. 

in particular, the invention has mainly been described in connection with 
GPS. However, it will be appreciated that tiie invention may also be applied 
20 to other signal reception systems where de-spreading of weak received 
signals is employed as a part of a signal detection system. For example, the 
invention may be applied to other positioning systems, e.g. the GLONASS 
system. 

25 Even though the invention has mainly been described in connection with a 
receiver in which the signal is mixed down to IF, it is understood that the 
invention may also be applied in connection with other receiver architectures. 
Examples of such architectures include an architecture in which the RF 
signal is directly down-converted to baseband, a so-called homodyne or 

30 zero-IF receiver. In another embodiment referred to as a superheterodyne 
receiver, the down-mixing is perfomried in several steps. 
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It is further understood that the signal acquisition described herein may 
advantageously be used in an autonomous GPS receiver. However the 
signal detection may also be used in a signal receiver receiving a priori 
information, e.g. a GPS receiver integrated in a mobile terminal or any other 
5 portable communications device where the mobile terminal receives 
assistance information from the mobile ne1wori<, e.g. information aiding the 
satellite selection, the frequency synchronisation, Doppler estimation, timing 
synchronisation, navigation message, or any subset or combination of the 
above. For example, in a network assisted system the method described 
10 herein may be used in situation without network coverage, to improve the 
indoor operation capability, and/or the like. In some embodiments in 
connection with a network-based system, the pseudo-range reporting to the 
network may be improved. 

15 The invention can be implemented by means of hardware comprising several 
distinct elements, and by means of a suitably programmed microprocessor. 
In the device claims enumerating several means, several of these means can 
be embodied by one and the same item of hardware, e.g. a suitably 
programmed microprocessor, one or more digital signal processor, one or 

20 more ASIC circuit, or a combination of the above. The mere fact that certain 
measures are recited in mutually different dependent claims or described in 
different embodiments does not Indicate that a combination of these 
measures cannot be used to advantage. 

25 It should be emphasized that the term "comprises/comprising" when used in 
this specification is taken to specify the presence of stated features, integers, 
steps or components but does not preclude the presence or addition of one 
or more other features, Integers, steps, components or groups thereof. 
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